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Abstract 

We explore the possibility of inducing in heterostructures driven by an ac gate 
voltage the coherent current suppression recently found for nanoscale conductors in 
oscillating fields. The destruction of current is fairly independent of the transport 
voltage, but can be controlled by the driving amplitude and frequency. Within a 
tight-binding approximation, we obtain analytical results for the average current 
in the presence of driving. These results are compared against an exact numerical 
treatment based on a transfer-matrix approach. 
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1 Introduction and modelling 

The study of electron transfer comprises a rich variety of systems in many 
different areas such as chemistry, biology, and life sciences [1,2]. Although 
electron transfer processes are mainly attributed to electrochemical applica- 
tions, they are conceptually related to molecular electronics [3-5] and electron 
transport in low dimensional materials in solid-state physics. In that context, 
semiconductor heterostructures represent a popular physical system for the 
investigation of mesoscopic transport [6-8] and tunnelling phenomena [9-13]. 
The main reason for this is the high mobility and the rather long mean free 
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path of the charge carriers populating them. Standard beam epitaxy tech- 
niques make the accurate growth of alloys of such materials on substrates 
possible, and the nearly identical lattice parameters, together with the possi- 
bility of controlling the band gap, turns the combination GaAs/Ala.Gai_a;As 
into an ideal candidate for building complex low dimensional structures with 
quantum wells and tunnel barriers. Moreover, these setups open various ways 
to study tunnelling in time-dependent systems [14-16]. A straightforward pos- 
sibility for introducing a time- dependence is the application of an ac transport 
voltage which only modulates the energies of the electrons in the leads while 
the potential inside the mesoscopic region remains time-independent. This 
kind of driving allows for a description within Tien-Gordon theory [17] which 
expresses the dc current in terms of the static transmission and an effective 
distribution function for the lead electrons. If the time-dependence enters via 
an external microwave field or an ac gate voltage, however, such an approach 
is generally insufficient [18]. 

A remarkable difference with respect to the static situation is the emergence 
of inelastic transport channels stemming from the emission or absorption of 
quanta of the driving field. For a periodically time-dependent transport situa- 
tion, however, we expect the transmission probabilities and, consequently, the 
resulting current to be time-dependent as well. This follows indeed from a re- 
cently presented Floquet theory for the transport through driven tight-binding 
systems [16, 18]. For the computation of the dc current, this approach justifies 
the applicability of a Landauer-like current formula where the static trans- 
mission is replaced by the time-averaged transmission of the time-dependent 
system. 

The transmission of both the elastic and the inelastic transport channels can 
depend sensitively on the driving parameters; the contribution of certain chan- 
nels can even vanish. For the transport across two barriers which enclose an 
oscillating potential well, Wagner [19] showed that it is possible to suppress 
the contribution of individual inelastic scattering channels. The total current, 
however, is given by the sum over all channels, and thus it is not possible 
to isolate the contribution of a single channel in a current measurement. By 
contrast, in the case of transport through a two-level system with attached 
leads, it has been found that driving with a dipole field has directly observ- 
able consequences. There, the driving not only affects the contribution of in- 
dividual transport channels, but the dc current can be suppressed almost 
entirely [20,21]. Therefore, for the appearance of this coherent current sup- 
pression, it is essential that the central region consists of at least two weakly 
coupled wells which oscillate relative to each other [18]. 

In this work, we explore the possibility of coherent current suppression in 
double-well heterostructures. Thereby, we compare two theoretical approaches 
to describe coherent transport in quantum- well structures: The transfer-matrix 
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Fig. 1. Model potential for the double- well heterostructure. In the numerical calcu- 
lations, we employ barriers with the heights = Vr = 90meV, Vq = 40meV and 
the widths c^l = c?r = 5nm, dc = 15 nm. The dotted lines mark the energy of a 
metastable tunnel doublet with splitting energy 2A. The left well is subject to an 
electric dipole field generated by an alternating gate voltage with amplitude 14c. 

method and a tight-binding approach. As a model we consider the triple- 
barrier structure sketched in Fig. 1 where the driving enters via an oscillating 
gate voltage which modulates the bottom of the left well. The applied trans- 
port voltage is assumed to shift the Fermi energy of the left lead by —eV with 
— e being the electron charge. We note that since the time-dependent gate 
voltage affects only one well, the structure depicted in Fig. 1 is sufficiently 
asymmetric to also act as an electron pump, i.e. to induce a non-zero current 
for eV — [16]. In this work, however, we focus on the transport properties 
in the presence of a finite bias voltage. 

For the exact numerical computation of the transmission probabilities, we em- 
ploy the transfer-matrix method developed by Wagner [22] , which is reviewed 
in Section 2. In Section 3, we introduce the related tight-binding system for 
which the transport properties can be calculated analytically within a high- 
frequency approximation scheme [21]. The predictions from the perturbative 
approach are compared to the exact solution in Section 4. 



2 Transfer-matrix method 



Following Landauer [23], we consider the coherent mesoscopic transport as a 
quantum mechanical scattering process. The central idea of this approach is 
the assumption that sufficiently far from the scattering region, the electronic 
single-particle states are plane waves and that their occupation probability 
is given by the Fermi function with the chemical potential depending on the 
applied voltage. The unitarity of evolution under coherent ac driving allows 
us to write the resulting currents as [24] 

I=lJdE [Tr^{E)U{E) - T^r{E)UE)], (1) 
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where T^(E) denotes the total transmission probabihty — i.e. sTimmed over 
transverse modes and outgoing inelastic channels — of an electron with energy 
E from the left lead to the right lead while Tlr(£^) describes the respective 
scattering from the right to the left lead. For time-independent conductors, the 
time-reversal symmetry of the quantum mechanical scattering process together 
with the energy conservation ensures Tr,l(£') = Tlr(£') such that, in the 
absence of a transport voltage, the current vanishes. This is not the case for 
a general time-dependent structure [16]. 

When the total Hamiltonian is time-periodic due to an external driving field, 
H{x,t) = H{x,t + T), one can apply Floquet theory [14,25,26]. It states that 
the corresponding time-dependent Schrodinger equation has a complete set of 
solutions of the form 

ijja{x,t) ^ exp{-ieat/h)ua{x,t), (2) 

where Ua{x,t) = Ua{x,t + T) denotes the so-called Floquet states, and the 
so-called quasienergies in analogy to the quasimomenta of Bloch theory. 

Owing to their time-periodicity, we can decompose the Floquet states into a 
Fourier series 

oo 

u^ix,t)= J2 e-'"™Ma,n(a;). (3) 

n=— oo 

The form (3) of the Floquet states suggests that during the scattering, an 
electron with initial energy E evolves into a coherent superposition of states 
with energies E' = E + nMl. The arbitrary integer n is referred to as the 
sideband index; the inelastic channels are called sidebands. We emphasise 
that, despite the existence of a band bottom, the summation over n in Eq. (3) 
is unrestricted. 

A proper calculation of the dc current through a time-dependent scatterer 
must now include these inelastic channels, and the Landauer formula of Eq. (1) 
has to be conveniently generalised to take them into account. Since the Flo- 
quet scattering states can be thought of as having been created from the 
orthogonal dc states by adiabatically switching on the driving, they too must 
be orthogonal, so that it is sufficient to sum over channels incident from both 
leads [24]. 

In the transfer-matrix method described below, the Floquet states are decom- 
posed into plane waves throughout the driven heterostructure, so that they 
can be appropriately matched with the scattering channels in the leads. This 
decomposition allows us to separate the time- and the space-dependent parts 
of the wave function, thus obtaining directly the time-independent probabili- 
ties that go into Eq. (1). 

For a spatially constant potential with a time-dependent gate voltage 14c(^)) 
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the solution of the Schrodinger equation is readily obtained to read 

+00 . J . 

xlj{E,z,t)^ Mz)exp\--{E + nhn)t-i(t){t)\ (4) 

n=-oo In ) 

with the accumulated phase 

cf>it)^U'dt'V,,it')^<t>{t + T). (5) 
n Jo 

Its time-periodicity follows from the zero time-average of the gate voltage. 

Neighbouring layers of a heterostructure may have different ac voltages applied 
in addition to different band-edges. As a consequence, the wave functions in 
Eq. (4) which solve the Schrodinger equation in each layer do not coincide in 
the general case. The solution for the complete system has to be constructed 
by matching the corresponding wave functions at the interfaces between layers. 
With this goal in mind, we assume that the wave function (4) is a solution of 
the Schrodinger equation for the time-dependent Hamiltonian 



H{z,t) = Hoiz) + eV,,{t) 

—Y^z^)^z+^^^^+'^-''^'''^ 

and that, moreover, ipn{z) is an eigcnfunction of the time-independent Hamil- 
tonian Hq with the spatially piecewise constant effective mass ■m{z), and has 
the general form 

V'n(-z) = Anex-p{knz) + S„ exp(-A;„2;) . (7) 

The wave vector 

= [2m{V - E - nhn)]^^^ (8) 

describes travelling as well as decaying waves (bound states) for complex and 
real values of respectively. The matching conditions at an interface follow 
from the fact that both the wave function and the flux have to be continuous, 
i.e. at 2; = ^0 



lim ip{z,t)= lim iplzjt) 

Z^Zq Z—*Zq 

Id Id 
z-^z+ m{z) dz z^zo m{z) dz ^ ' 

This yields an infinite system of algebraic equations for the coefficients An and 
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in each layer. Inserting the Fourier expansion of the phase in Eq. (5), 

exp|--^V;cSinQt| = ^ J„/ (^^^j exp(-in'Qi), (10) 



where J„' is the n'-th order Bessel function of the first kind, allows one to 
recast these equations for an interface between layers I and II at 2; = Zi in 
matrix form: 



T 



Zi 



41 \ / 4" 



Bl '\B, 



The matrices T^. and with elements Tl..^^^, and T";„_„/, respectively, are 
of infinite dimension, and contain the coefficients of all possible scattering 
channels, i.e. photon exchanges between the incoming electron and the driving 
field, at either side of the interface. Their precise form depends on whether or 
not a layer is affected by the time- dependent gate voltage. The transfer matrix 
T^.^x between two sides of a layer of width then defined as [27] 



T 



^3 



^UTy-\ (12) 



This definition in terms of layers, rather than using a similar one across inter- 
faces, is physically more sensible, since the matrices depend on the properties 
of just one layer. They are also numerically easier to implement, as there are 
fewer qualitatively different matrices to deal with. To calculate the total trans- 
fer matrix across the structure, we have to multiply the matrices across all 
different layers and obtain 

Tl->r = T^j^T^^.^^j^ • • • T^L^z^T^L, (13) 

where T^^^ and T^j^ represent the initial and final matrices at the ends of 
the hctcrostructure. With the elements T^fl^i ^hc total transfer matrix we 
can find the probability that an electron with energy E + nhfl in lead L is 
scattered into a channel with energy E' — E+n'fUD, in lead R, with integer n, n'. 
The diagonal elements T^l!^^ are closely related to the (static) transmission 
probability Trl, while the off-diagonal elements T^^^ describe the effects of 
the absorption or emission oi n — n' photons on the transmission probability 
of the electron. For fiat conduction bands on both sides of the heterostructure, 
the wave functions in the contacts are plane waves, and in this case the proper 
boundary conditions to describe an electron incident from, say, the left-hand 
side at energy E are — 5„^o and — 0. The transmission probabihty in 
sideband n is then defined as 

_ A£mL A| ^ 
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Fig. 2. Schematic energy diagram of a symmetric double-well structure in the 
tight-binding approximation. The tunnelling matrix element is given by A and each 
well couples with Fl and Fr to the associated lead. The state |L) in the left well 
is driven by a oscillating gate voltage with amplitude T^c In addition, an external 
bias V = (/xr — /UL)/e is applied. 

where k^^ and represent the wave vectors on the right- and left-hand side 
in sidebands n and 0, respectively. 



In a numerical implementation of the transfer-matrix technique, it is necessary 
to truncate the infinite matrices. Thereby for consistency, a proper cut-off has 
to be so large that unitarity of the scattering process is preserved, i.e. 

-t-oo -foo 

E E ^LL = 1, (15) 

n=— oo n=—oo 

where i?LL represents the reflection probability of an electron to be reflected 
from energy E into a sideband at energy E + nhfl on the same side. The 
number of sidebands that need to be taken into account to meet a given 
accuracy (which in our calculations was set to 10~^^) depends essentially on 
the ratio eVg^/Ml, as this is the argument of the Bessel function that 
determines the weight of each sideband n. To proceed, one starts at the initial 
value of \4c = with a tentative number of sidebands, and increases it for 
growing driving amplitudes if a check with Eq. (15) suggests that unitarity is 
breaking down. When particle number conservation is restored one can go to a 
higher T4c- Transfer matrices such as those employed here have the advantage 
of being easily scalable to arbitrarily complex structures. The combination of 
flexibility in structural properties and numerical accuracy makes this method 
well-suited to the study of strongly driven semiconductor heterostructures. 



3 Tight-binding approximation 



A different approach to study resonant tunnelling in a driven double-well 
structure is based on the adoption of a tight-binding (TB) picture. Figure 2 
depicts the TB configuration corresponding to the heterostructure introduced 
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above. 



3.1 The model 

Then, the Hamiltonian of the system is given by 

H{t) — H^ellsit) + -f^leads + -f^contacts , (16) 

where the individual terms describe the driven quantum wells, the electron 
reservoirs in the leads, and the coupling of the left and the right well to the 
respective neighbouring lead. For simplicity, we neglect the electron spin. 

Within the framework of the TB approximation, the time-dependent quantum- 
well Hamiltonian reads 

-f^weiis(i) = X] Hi^e'{t)clce' = -A(4cr + ^cl) + eKc cos(Qi)c[cL . (17) 

An electron can be localised in the left or right well, whereupon the fermion op- 
erator Q (cj) annihilates (creates) an electron in the respective well {i = L, R). 
These localised states |L) and |R) are coupled by the tunnelling matrix ele- 
ment A. For convenience, the energy scale is set such that the on-site energies 
of the two resonant TB levels are zero and lie exactly halfway in between 
the transport bias window defined by the chemical potentials //l and /xr. The 
second term of the Hamiltonian (17) accounts for the harmonic driving of 
the traversing electrons in the left well via an oscillating gate voltage with 
amplitude Vg^ and period T — 2Tr /fl. 

The leads are modelled as ideal Fermi gases with the Hamiltonian 

Hleads = ^^iqclqCiq, (18) 

e,q 

where C£q (cl^) annihilates (creates) an electron in the lead with energy e^g with 
i = L,K. As an initial condition, we employ the grand-canonical ensembles of 
electrons in the leads at inverse temperature (3 = l/Zc^T. Therefore, the lead 
electrons are characterised by the equilibrium Fermi distribution ft{eiq) ~ 
{l + exp[-/3(e^,-/x^)]}-i. 

The localised state in each well couples via the tunnelling matrix element Vnq 
to the state \iq) in the respective lead. The Hamiltonian which describes this 
interaction has the form 

-f^contacts = XI ^(■lAq'^t + ^-C. (19) 
t,q 
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The lead-well coupling is entirely specified by the spectral density r^(e) = 
27r |V£qp5(e — efg). Since, for the system at hand, the bandwidth of the 
conduction band of the leads is much larger than the energy regime where 
transport happens, the spectral densities are practically constant, i.e. T^i^e) ~ 
r^, which defines the so-called wide-band hmit. 



3.2 Floquet transport theory 



In order to determine the time-averaged dc current which matches Eq. (1) for 
the TB approximation, we employ a generalised Floquet approach to solve the 
corresponding Heisenberg equations of motion and derive an expression for the 
retarded Green's function in terms of Floquet states. This result is apphed to 
study transport by evaluating the operator /^(t) = {ie/h)[H(t),J2qclgCiq] for 
the time-dependent current through contact i. In the next subsection we will 
also derive a analytic expression for the current valid in the high frequency 
limit with respect to the driving. 

We start out by stating the solution of the Heisenberg equations of motion for 
the lead operators, which is given for the left lead by [18] 

CLq{t) = CL,(to)e-'^^^(*-*°)/'^ - '^'jAt'e-''-^^'-''^'^c^{t'). (20) 

to 

For the corresponding solution of CRg(i), L has to be substituted by R. Now 

inserting these solutions into the Heisenberg equations of motion for the quan- 
tum wells yields the two coupled linear differential equations 

ic r i 

CL = - ^Kc COs(l]t)cL - ^CL + ^ACR + ii^{t) (21) 

i r 

CR=-^AcL-^CR + eR(i). (22) 

Here, within the wide-band limit, the coupling to the leads has been eliminated 
in favour of the spectral density Ti and the fermionic fiuctuation operator 

Ut) = E ^4 e-^^''(*-*°)/'^Q,(io). (23) 
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Assuming that the initial conditions are those of the grand canonical ensemble, 
this Gaussian noise operator satisfies 

Ut)) = 0, (24) 
iClit) Ce{t')) = S,e'^, Jdee'<'-''y^fe{e). (25) 
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Then the operator for the time-dependent current through the left lead be- 
comes [18] 



m = ^c[(t)cL(t) - e[ci{t)Ut) + d{t)cL{t)] (26) 

with a corresponding expression for Iji{t). Here we made use of the wide- 
band limit. To evaluate the time-dependent current, we thus have to find the 
sohition for the inhomogeneous set of the quantum Langevin Eqs. (21) and 
(22) of the quantum-well operators, which is formally given by 

ce{t) = j^dT [G',L(t, t - T)i^{t - r) + G,^{t, t - T)U{t - r)] (27) 

in the stationary limit to ^ oo with i — L,R. What remains is to determine 
the retarded Green's function G{t, t — r). This is where Floquet theory comes 
into play by making use of the T-periodicity of the driving. Solving the Floquet 
eigenvalue equation 

n^eiisit) - iS - ih^^ \ua{t)) = (e„ - ihja)\u^{t)) (28) 

of the physical problem at hand, where H{t) — J2e,e' \'^)He,e'{t){i'\: we get 
the Floquet states \ua{t)) and the complex- vahied quasienergies — ihja- 
Note that the prior equation is, in contrast to the usual Floquet equation, 
non-Hermitian. This is due to the presence of the self-energy 2S = |L)rL(L| + 
|R)rR(R|, which results from tracing out the leads [28]. Therefore, Eq. (28) has 
to be solved also for its adjoint eigenstates {w^it)) [16]. With the corresponding 
expression for the propagator U {t, t—r), the retarded Green's function assumes 
the form 

G{t,t- t) = Ee-'^^"/'-'""^1«a(t))«(t - r)|e(r), (29) 

a 

where ©(r) is the Heaviside step function. 

The dc current is now obtained by calculating the expectation value (/l(^)) 
and averaging over one driving period. This time-averaging will cancel those 
terms of {Ii,{t)), which are responsible for a T-periodic charging of the wire. 
After eliminating backscattering terms [18], we arrive at the very compact 
form for the final result 



e 
h 



E lTS{e)Me) - T^i{e)Me)], (30) 

n •' 



where the total transmission probabilities are given by (e) = rLrR|GL'^(e) p 
and T^{e) = TiXR\G^^^{e)\'^, which resemble the Fisher-Lee relation [29]. The 
above equation also holds for the current through the right contact owing to 
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charge conservation. The Green's function 



^''-i;:,e-ie^ + n'hn-ih^^) ^^^^ 

is the t-averaged Fourier transformed of the propagator (29). PhysicaUy, it 
describes the propagation of a transmitted electron with initial energy e from 
one lead to the other lead undergoing scattering events with emission (n < 0) 
or absorption {n > 0) of |n| photons, or being transmitted elastically {n — 0). 



3.3 High-frequency limit 



The Floquet treatment of the present transport problem allows for the imple- 
mentation of a stationary perturbation scheme for driving frequencies much 
larger than all other frequency scales of the system [26]. This approach has 
recently been extended to transport situations which are characterised by 
the presence of leads [18, 21]; here we only outline the derivation and re- 
fer the reader to Ref. [16]. A particular benefit of this perturbation scheme 
is the achievement of a physical understanding of the transport processes 
by approximately mapping the time-dependent problem to a static one with 
renormalised parameters. For the static situation, in turn, the current is well 
known by Eq. (1), where the transmission in the wide- band limit reads T(e) = 
^LR(e) = rLrR|G'LR(e)p. This looks similar to the driven case but with n = 
and therefore we have in contrast to the driven system GLR(e) = GRL(e)- For 
the static system given by (16) setting Vj^c = 0, we obtain for the transmission 

assuming equal coupling to the leads (F^ = F). 

The starting point of the approximation scheme is the unitary transformation 
Uo{t) = exp |--^V;cSin(Qi)ctcL| , (33) 

which is first applied to the quantum- well Hamiltonian (17). For sufficiently 
large driving frequencies Q ^ A/h, a separation of time scales is performed by 
this transformation. Thereby, fast oscillations of the transformed Hamiltonian 
are neglected by averaging over a driving period [14,30]. Finally, we arrive at 
the effective Hamiltonian for the quantum wells 

^eff = ^ £dt (c/o^ H^,i^{t) Uo - ihUlUo) 

= -Aeff(c[cR + 4cL), (34) 
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which is of the same form as in the static case but with the effective tunnelhng 
matrix element A^s = ^o(eKic/^^)^ [20,30]. Jo is the zcroth order Bessel 
function of the first kind. Therefore, the resulting effective transmission Teff(e) 
with the substitution A — > Aeff in Eq. (32) is controllable via the driving 
parameters and even vanishes at zeros of the Bessel function. 

The transformation (33) also affects the lead-well coupling. If we apply Uo{t) 
also to i^contacts and solve the Heisenberg equations for the lead and quantum- 
well operators in the wide-band limit, we can eventually extract the new fluc- 
tuation operator. For the left lead one finds 



whereas ijRlt) remains unaffected. Now calculating the correlation function 
and time-averaging it over one driving period to neglect the T-periodic con- 
tributions, the resulting expression assumes the form (25) but with the Fermi 
function of the left lead replaced by the effective electron distribution 



The squares of the nth-order Bessel function of the first kind in this expression 
weight those processes where an electron with energy e is transmitted from the 
left lead to the double-well system under the emission (n < 0) or absorption 
{n > 0) of \n\ photons. The effective electron distribution exhibits steps at the 
energies e = /^l + nhfl and is constant elsewhere. 

With the effective quantities (e) and /L,efr(e) the driven problem is ascribed 
for fast driving to a static one. Since Tcff(e) is sharply peaked around e = 
and /r(0) and /L.cff(O) are constant for finite voltage, the current in the high- 
frequency approximation results in 



applying the effective parameters to the current formula (1). Here K(y) is a 
shorthand notation for the integer part of e\V\/2hfl. 

In order to compare the transfer-matrix and the tight-binding approach, we 
have to ensure that the same physical situation is addressed. As a matching 

condition we compare the transmission T(e) in the time-independent case 
(K,c = 0). The level splitting energy 2 A due to the central tunnel barrier is 
extracted from the resonance peaks of the doublet states computed within the 
transfer-matrix method. Solving for F in Eq. (32) with e = and T(0) and A 





(36) 




(37) 
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Fig. 3. Average current vs. driving amplitude obtained numerically from trans- 
fer-matrix (solid line) and tight-binding (dashed) methods. Also shown is the 
high-frequency approximation (dashed-dotted) . The inset depicts the value of the 
first current minimum as a function of the driving frequency. Solid (transfer-matrix) 
and dashed (tight-binding) line decay approx. as 1/0. The chosen parameters are 
m = 1.15 meV, V = 6.0 mV, T = 0.16 meV and A = 0.23 meV. The corresponding 
parameters for the barriers are the same as those of Fig. 1. 

taken from the previous calculation, the corresponding lead-well coupling for 
the tight-binding system is determined. 



4 Coherent transport suppression 

We now turn our attention to the coherent control of current. Tunnelling sup- 
pression in a closed, driven system is known for more than a decade. For exam- 
ple for a driven bistable potential, tunnelling breaks down at exact crossings of 
the quasi-energy spectrum, i.e. one observes the so-called coherent destruction 
of tunnelling [31,32]. Tunnelling suppression has been studied in a number of 
cases [22,33,34], but the investigation in a transport context, i.e. in an open 
system where an appropriate treatment of the leads is crucial, has received 
attention only recently [18,21]. 

Surveying the time- averaged current calculated numerically from the transfer- 
matrix and the tight-binding method plotted in Fig. 3, we observe current 
minima for distinct values of e\4c/^^ fo^' frequencies in the microwave regime. 
The reason for the current suppressions becomes apparent by comparison with 
the high-frequency approximation, which exhibits minima close to those of the 
transfer-matrix and tight-binding curves. The current (37) vanishes whenever 
the ratio eVac/M), assumes a zero of the Bessel function Jq, i.e. for the values 
2.405, 5.520, 8.654, . . . , since then A^s oc Jg = 0. By varying the ratio between 
driving amplitude and frequency, we can thus tune the tunnelling between the 
two wells and thereby control the current [21]. For a frequency Q = 5A/h, 
the analytical expression (37) shows a remarkable agreement with the exact 
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Fig. 4. Deviation of the driving amplitude for the first current minimum from the 
expected first zero of Jq, Vq = 2.405 ^Jl, for different barrier widths and heights. 

The parameters for the first three data sets arc dc = 15 nm, p, = 12.0 meV and 
Vc = 40 meV, whereas for the last one (o), we chose p, = 13.3 meV and Vc = 80 meV. 

tight-binding result (30) for T4c ^ V. The inset of Fig. 3 shows the minimum 
current at the first suppression decays as a function of the driving frequency 
Q. This is expected from the good agreement between the numerical results 
and the high-frequency approximation, because the approximation accounts 
for the first order term in a perturbative scheme in 1/^2 [18,21]. Higher order 
contributions are included in a numerically exact calculation, which results in 
a non-vanishing current at the minima. A similar Q-dependence is observed 
also for the transfer-matrix formalism. 



While the general shape and magnitude of the current are very similar for both 
models, there still appears a small difference in the location of the minima for 

the relatively low barriers chosen in Fig. 3. For a continuous potential, the 
current assumes minima at values of eV^c/fi^ higher than those predicted by 
the tight-binding description. We can understand this shift by analysing, for 
given fl, the deviation 5V — Vmin — Vo of the driving amplitude Vuiin for which 
the current exhibits its first minimum. The amplitude Vq corresponds to the 
first zero of Jq. In Fig. 4 we plot SV as a function of nd, where d = di^ = d^i 
and K = [2m{V — fL)/h^]^^'^, i.e. nd is the instanton action in units of h and 
exp{—2Kd) is the WKB transmission probability of the outer barriers in Fig. 1. 
Here V — Vi, — is the corresponding barrier height and jj, — {fxi, + fXR)/2 
denotes the average chemical potential representing approximately the mean 
energy of the resonance doublet. 

If the width of the outer barriers is kept fixed, 6V decreases for growing nd 
because then the resonance energies are further away from the barrier edge. 
Therefore, the wave functions of the well states become more localised. This 
situation corresponds in the tight-binding picture to a lead-well coupling F 
that is almost energy independent and thus reproduces the wide-band limit. 
Furthermore, this argument is used to explain the smaller deviation observed 
with thinner barriers for the same nd, since Vl is much larger in that case. 
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As can be seen by comparing data sets for different central barrier heights in 
Fig. 4, an increase of the height of the central barrier Vc reduces the level split- 
ting 2A, that is, the overlap between the localised states in the left and right 
well in a tight-binding description. Thus the tight-binding and the transfer- 
matrix results converge as a function of the barrier height. Finally, it is impor- 
tant to note that varying any of the barriers affects the transmission properties 
of the whole heterostructure, in contrast to the tight-binding model, where the 
different coupling parameters can be controlled independently. 



5 Conclusions 



We have demonstrated that the dc current across a double-well heterostruc- 
ture can be suppressed by the purely coherent influence of an oscillating gate 
voltage. We have used a transfer-matrix method as an exact approach to com- 
pute tunnelling currents through such a system. We compared these results to 
that obtained from a tight-binding Floquet description. In particular, we find 
that the current suppression is controlled by the ratio of the driving frequency 
and amplitude. This can be understood by exploring the high-frequency limit 
within the tight-binding formalism. In this perturbative scheme, the time- 
dependent system is mapped onto a static one with renormalised parameters, 
that is, with an effective hopping matrix clement accounting for inter-well 
coupling and with an effective electron distribution for the attached left lead. 
Since the effective inter-well coupling depends on the ratio between driving 
amplitude and frequency, the transport properties of the double well can be 
adjusted through the driving parameters, with the effective behaviour ranging 
from transport through an almost open channel to a regime of rare tunnel 
events. 



The results presented in this work strongly support the idea that transfer- 
matrix and tight-binding descriptions of quantum transport are equivalent 
provided the barriers are sufficiently high. In this case the lowest resonance 
states in the wells are rather localised and consequently the tight-binding de- 
scription becomes accurate. For a proper choice of parameters, we find a good 
agreement between the exact transfer-matrix calculation and the results ob- 
tained within the tight-binding formalism. The study presented here shows 
that the coherent control of time- dependent electron transport can be inves- 
tigated with current semiconductor nanotechnology. 
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